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Abstract: This work aims to present the current state of simulations of electroluminescence (EL) 
produced in gas-based detectors with special interest for NEXT — Neutrino Experiment with a 
Xenon TPC. NEXT is a neutrinoless double beta decay experiment, thus needs outstanding energy 
resolution which can be achieved by using electroluminescence. The process of light production 
is reviewed and properties such as EL yield and associated fluctuations, excitation and electrolu- 
minescence efficiencies, and energy resolution, are calculated. An EL production region with a 
5 mm width gap between two infinite parallel planes is considered, where a uniform electric field is 
produced. The pressure and temperature considered are 10 bar and 293 K, respectively. The results 
show that, even for low values of VUV photon detection efficiency, good energy resolution can be 
achieved: below 0.4 % (FWHM) at Q^^ = 2.458 MeV. 

Keywords: Scintillators, scintillation and light emission processes (solid, gas and liquid 
scintillators); Large detector systems for particle and astroparticle physics; Time projection 
chambers; Detector modelling and simulations II (electric fields, charge transport, multiplication 
and induction, pulse formation, electron emission, etc). 
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1. Introduction 

The Neutrino Experiment with a Xenon TPC (NEXT) will search for the j8j8°^ decay in ^^^Xe using 
a 100 kg, high-pressure gaseous xenon, electrominescent time projection chamber (HPXe TPC). 
The project is approved for operation in the Canfranc Underground Laboratory (LSC), Spain. The 
TPC will have separated readout systems for calorimetry and tracking to facilitate both measure- 
ments (see Figure |I]). 

The two electrons emitted by the jSjS"^ decay transfer their energy to the medium through ion- 
ization and excitation with the characteristic topological signature "spaghetti with two meat balls". 
The excitation energy results in the prompt emission of vacuum ultraviolet (VUV) light - primary 
scintillation. The ionization tracks (positive ions and free electrons) left behind by the particles are 
prevented from recombination by a suitable drift electric field. Negative charge carriers - primary 
electrons - drift toward the TPC anode, entering a region with a more intense electric field created 
by two parallel meshes (EL region in Figure There, they are accelerated and collide with the 
gas atoms. Depending on the intensity of the electric field, electrons can either elastically collide 
with the atoms, or excite or even ionize them. As a result of the de-excitation of excited atoms, 
further VUV photons are generated isotropically - secondary scintillation or electroluminescence 
(EL), Section Therefore, in such an electroluminescent TPC, both the scintillation and ionization 
processes ultimately produce VUV photons, to be detected with a photosensors array (photomul- 
tiplier tubes, PMTs) located behind the cathode. The detection of the primary scintillation light 
constitutes the start-of-event, to, whereas the detection of EL provides an energy measurement. 
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Figure 1: Schematic representation of the Neutrino Experiment with a Xenon TPC detector and of 
the topological signature "spaghetti with two meat balls" of the two electrons emitted by the jSjS"^ 
decay. 



Electroluminescent light is used also for tracking, as it is detected also by a second array of pho- 
tosensors (e.g. multi-pixel photon counters, MPPCs) located behind and close to the EL region 

t- 

Electroluminescence has already demonstrated in the past to be a very good option in order to 
achieve excellent energy resolution [§, ||], which is of particular importance in neutrinoless double 
beta decay searches. As an important example, high resolution spectrometry of low-energy X-rays 
is possible with Gas Proportional Scintillation Counters [0, ||]. 

Because of this, it is important to understand the process of secondary light emission (electro- 
luminescence) and try to assess its different properties, namely excitation and electroluminescence 
efficiencies, EL yield and corresponding fluctuations. With these parameters in hand, it is pos- 
sible to predict the light gain and energy resolution that are achievable with a particular detector 
geometry. This assessment can be done through detailed simulations as explained in reference [^]. 
In this work we present results obtained applying detailed Monte Carlo simulations to the NEXT 
experiment conditions. 

2. Electroluminescence process 

At pressures above few hundred Torr [Q, the main channel of de-population of the excited atoms 
is through the fonnation of excimers — electronically excited molecular states. Excimers, are 
formed through three-body collisions between one excited atom, R*, and two atoms in the ground 
state, R: 

R* + 2R^R*2* +R (2.1) 
The excimers are formed through process ( [2. 1| ) in high vibrational states and, at the pressures 
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studied in this work, collide with ground atoms, R, losing vibrational energy: 

R*2*+R^R*2 + R (2.2) 

The resultant excimers, /?2, in a low vibrational state, emit VUV photons: 

R*2^2R + hv (2.3) 

This results in the so-called "second continuum" centered at 173 nm and with a FuU-Width- 
at-Half-Maximum (FWHM) of 14 nm Further details about the electroluminescence process 
in pure noble gases can be found in Ref. [^. 



3. Simulation 



3.1 Toolkit 



The platform used to perform the simulations is the new C++ version of the microscopic technique 
of GARFIELD [^]. This is a Monte Carlo technique which tracks the electrons at the atomic level 
using, currently, procedures and cross sections available in MAGBOLTZ 8.9.3 [|l^, 11]. 



The electron follows a vacuum trajectory between collisions, the path length being sampled 
according to an exponential distribution around the (energy-dependent) electron mean free path. 
The null-collision technique is used to consider the variations of the mean free path due to changes 
of the electron kinetic energy between collisions [|^]. In pure xenon, the program classifies the 
collisions as elastic, inelastic (electronic excitation) or ionization. It is possible to obtain informa- 
tion about each excited atom produced in the gas, namely: spatial position, time, excitation level 
and energy spent in the excitation. 

Figure ^ shows the cross sections used in our simulations. Detailed simulation of xenon exci- 
tations is possible by parametrizing the excited energy levels as a function of 50 energy groups. 

An anisotropic angular distribution for elastic collisions is implemented by using both the 
momentum transfer (EMT) cross section and the total elastic (ELT) cross section according to the 



method described in Ref. []13[]. 
3.2 Model 

Each primary electron is allowed to drift a distance of d = 5 mm under the influence of a uniform 
electric field. The effect of the field distortion near the meshes is not considered since it has no 



significant effect []15|]. We assume that the gas is at a pressure of p = 10 bar and at a temperature of 
293 K. A set of A'^e = 40, 000 primary electrons is used for each value of the voltage applied between 
the parallel planes, V. The starting direction of each primary electron is generated isotropically. 
The starting energy is generated according to the energy distribution calculated by Magboltz for 
the actual reduced electric field, (^^j = j^, (see Figure ^. 

We assume that every excited atom gives rise to the isotropic emission of a VUV photon. Its 
energy is generated according to a Gaussian distribution with the characteristics of the "second 
continuum" (see Section 

The model was previously validated in Ref. ^ by comparing the simulations with available 
values of the EL yield, the excitation and the electroluminescence efficiencies. Earlier Monte Carlo 
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Figure 2: Cross sections used by Magboltz 8.9.3 for Xe. The program includes detailed cross 
sections for elastic collisions, excitations (represented in 50 energy groups) and ionizations. 
Anisotropic elastic scattering is considered by using both the momentum transfer (EMT) and the 
total elastic (ELT) cross sections, shown as continuous and dashed gray lines, respectively. The 
blue lines refer to the first four excited states of Xe corresponding to the electronic configuration 



5p^6s' (Racah notation is used [14]), the green lines correspond to the 10 levels of the 5p^6p^ con- 
figuration, and the gray dotted lines to other sparsely populated levels. The red line corresponds to 
ionization (ION). 



results and experimental measurements [17| in uniform electric field geometry are presented 
in Figure ^ to illustrate the good agreement. 



4. Results 



4.1 Efficiencies 



Excitation and electroluminescence efficiencies, Qexc and Qel, respectively, as a function of the 
reduced electric field are shown in Figure |[ The first quantity represents the fraction of energy that 



is supplied by the electric field to the primary electrons and that is spent in excitations Jig ] : 



Q = tl (4.1) 

edKE 



In Eq. (|4. ID, ?iexc = 50 is the number of excitation groups available in Magboltz 8.9.3 for xenon, 
n' is the number of excitations of group / produced by the set of Ng primary electrons, gg^c the 
energy of the excitation group, E is the modulus of the applied electric field and e is the charge 
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Figure 3: Energy distributions of the electrons before each collision with the xenon atoms for 
different values of the reduced electric field. I 
broader and the peak shifts to higher energies. 



different values of the reduced electric field. As the value of f ^ j increases, the distribution gets 



of one electron. The second parameter, Qel, is the ratio between the energy emitted in the form of 
VUV photons and the energy transferred from the field to the A'e electrons: 



Gel 



edNeE 



(4.2) 



In Eq. (p^), is the energy of the VUV photon emitted after the excited atom of the group 
is produced. The VUV photon energy is randomly generated according to Section ^parameters. 

We note that 2el is always lower than 2exc- This is due mainly to the loss of vibrational 
energy by excimers before they emit a VUV photon (as explained in Section ^). The fast increase 
of both efficiencies for > 5.0 kVcm^^bar^^ is due to the contribution of secondary charges 
that start to be produced at high reduced field values. 



4.2 Electroluminescence yield 

The reduced electroluminescence yield, as a function of the reduced electric field is also 

shown in Figure ^ The reduced electroluminescence yield is defined as the average number of 
VUV photons emitted per primary electron and per unit of drift length divided by the pressure 
of the gas, p. The behaviour of with is approximately linear even if some ionization 
occurs. This happens for a < 0.1 ions • cm^\ being a the First Townsend coefficient (see Figure 
defined as the number of secondary charges produced when one electron crosses a distance of 
1 cm. 
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Figure 4: Reduced EL yield, Y/p, as a function of the reduced electric field, E/p. Excitation 
efficiency, Qexc, and EL efficiency, Qel, are also shown. Full symbols are the results of this work. 
Earlier Monte Carlo results for 2exc and Qel Jl^], as well as experimental measurements of the 



reduced EL yield [17], are also shown (open symbols) for comparison. 



Performing a linear fit to the simulation results, we obtain the following dependence: 
= (136 ± 1) - (99 ±4) [photons electron"' cm"' bar"'] 



(4.3) 



where E /p is expressed in kVcm 'bar '. Eq. (4.3) is in good agreement with experimental data 
measured at 1 bar [|T7]|: 



1401 



(4.4) 



^ — 116 [photons electron ' cm ' bar '] 
4.3 Fluctuations and energy resolution 

The FWHM energy resolution. Re, of an EL detector (such as NEXT) corresponds to the quadrature 
sum of the contributions of the different processes happening in the detector: 



Re = 2V21n2, 



^ 1 



\nI N,\nIe, 



+ 



+ J- 1 = 



ep 



ep 



Go 



(4.5) 



In Eq. (4J) the factor l^Jllnl corresponds to the relation between the FWHM and the stan- 



dard deviation, a, of a given probability distribution (FWHM = 1^/1 \nlo ~ 2.35a). The first term 
of the expression is related to fluctuations in the number of primary charges created per event (that 
is a jSjS"^ decay in NEXT), N^, the second to fluctuations in the number of EL photons produced 
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Figure 5: Simulated first Townsend coefficient, a , as a function of tlie reduced electric field. 



per primary electron, A^el. the third reflects the variations in the number of photoelectrons extracted 
to the PMT photocathode per decay, N^^, and the fourth the distribution in the the PMT's single 
electron pulse height, Gq. 

The primary charge fluctuations are well known and described by the Fano factor, F = a^/Ne. 
In this work we consider f = 0.15 [19]. In literature are reported values between 0.13 and 0.17 for 
pure xenon [20, n), |2[ 24] which we used to estimate the overall energy resolution error bars 
presented in Figure |8} 

The fluctuations associated with the electroluminescence production are described by the pa- 
rameter J defined as the relative variance in the number of emitted VUV photons per primary 
electron: 



"EL 

Nel 



(4.6) 



Concerning the third term, the conversion of VUV photons into photoelectrons follow a Pois- 



son distribution and thus 



ep- 



The fluctuations in the photoelectron multiplication gain within the PMT can typically be 
described by ( ^ ) =1 [l25|]. Taking into account the previous relations, Eq. (031) can be rewritten 
as follows [pi : 



2 



(4.7) 



ep 



Usually J is much smaller than F and the energy resolution of a EL detector is attributed only 
to the fluctuations in the primary charge production and in the photon detection system (a plane of 



PMT's in the case of NEXT) - first and third terms of Eq. ( |4-.7[ ). Given that we deal with simulated 
data, in this work we can also calculate the parameter / for the present geometry. Figure S shows / 
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Figure 6: Relative variance in the number of emitted EL photons, J, as a function of the re- 
duced electric field. It is also shown the Fano factor - the relative variance in the number of 
primary charges produced per event - for high pressure xenon and a total deposited energy of 
Qpfi = 2.458 MeV. 



as a function of the reduced electric field. As the electric field increases, the parameter / decreases 
until secondary electrons begin to be produced and the avalanche fluctuations start to dominate. For 
J, values that are much smaller than the Fano factor, can be achieved in optimal conditions, namely 
for reduced electric fields between 1.5 and 3.5 kV cm^^ bar^^ In order to better understand the 
effect of avalanche fluctuations in EL fluctuations, in Figure we show the A^el distributions for 
different values of the reduced electric field. The effect of secondary charges is clearly visible, 
since a long tail in the distribution appears for higher values of the electric field (Figures ^ and 



7d) 



Using the data obtained by the simulation, we can estimate the energy resolution which can be 
achieved by the NEXT detector. We assume that all the energy of the two j8j3"^ decay electrons, 
Qfifi = 2.458 MeV, is deposited in the gas medium, with an average energy lost per ion pair formed 



of w = 21.9eV [19, 24]: 



- Qrr 2.458 X 10^ eV 

= =^ = ^ 1 12, 237 electrons (4.8) 

w 21.9 eV 



Taking into account the intensity of the electric fields in the drift and electroluminescence 
regions that are planned to be used, we assume that all the primary electrons arrive to the EL 
region and cross the first mesh without either any recombination or attachment to electronegative 
contaminants in the drift region. 

The average number of photoelectrons produced in the plane of PMTs per decay, Nep, can be 
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Figure 7: Distribution in the number of emitted EL photons per primary electron, A^el. during the 
drift along d = 5 mm of xenon at 10 bar and 293 K, for different values of the reduced electric 
field. 



obtained as: 

Fep = kNeNEL (4.9) 

where k is the fraction of EL photons produced per jSjS"^ decay that gives rise to the production of 
a photoelectron. 

The FWHM energy resolution, Re, is shown in Figure || as a function of the reduced electric 
field for two different scenarios: an optimistic (k = 0.06 [|T9|]) and a conservative (k = 0.005 [[T]]). 
The contributions of each of the three terms of Eq. ( p~7| ) are also shown. The optimal energy 
resolution is achieved for a value of that is higher than the one marking the beginning of 
ionization, shown in Figure ^ as the reduced field value minimizing the value of /. The reason 
is that the still small additional fluctuations introduced by the production of secondary charge are 
compensated by the increase in the number of EL photons produced per primary electron due to 
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Figure 8: FWHM energy resolution, Re, as a function of the reduced electric field for two different 
scenarios: an optimistic (circles, k = 0.06) and a conservative (squares, k = 0.005). The values of 
Re if only each of the three terms of Eq. ( p7| ) would contribute are also shown. 



the increase of the EL yield. Also, and as expected, the reduced electric field value which gives the 
best energy resolution increases as k decreases, and the best value of Re is worst for lower values 
of k. For the experimental conditions considered here, the best energy resolution is obtained for 
between ~ 4.0 and ~ 5.0 kVcm^'bar^', that is the value usually attributed to the ionization 
threshold of Xe [^. 

5. Conclusions 

The C++ version of the microscopic technique of Garfield, allows to gather information on the 
excited atoms produced during the drift of electrons through a xenon gas detector. Such simula- 
tions rely on detailed cross sections for elastic collisions, excitations and ionization as modelled by 
Magboltz 8.9.3. By considering that every excited atom leads to the emission of a vacuum ultravi- 
olet (VUV) photon, we are able to accurately simulate electroluminescent (EL) hght production in 
realistic detector geometry setups. 

In this work we consider EL light production in the NEXT j8j8°^ experiment. We evaluate the 
EL yield and the corresponding fluctuations produced in such a setup, corresponding to a EL region 
given by a 5 mm uniform electric field gap of xenon at 293 K temperature and 10 bar pressure. We 
obtain that a FWHM energy resolution of 0.4% or better at the Q value of the jSjS"^ reaction, Qpp = 
2.458 MeV, can in principle be obtained, even considering conservative assumptions regarding the 
VUV photon detection efficiency. 

In conclusion, our simulations indicate that the 1% FWHM energy resolution goal of the 
NEXT experiment can in principle be met and even be improved upon, as far as fluctuations in 
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the primary ionization charge, in the electroluminescence, and in the VUV photodetection pro- 
cesses are concerned. This energy resolution goal is one of the main figures of merit for the NEXT 
experiment 
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